Linear shear flow past a hemispherical droplet adhering to a solid surface 
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This paper investigates the properties of a three dimensional shear flow overpassing a hemispher- 
ical droplet resting on a plane wall. The exact solution is computed as a function of the viscosity 
ratio between the droplet and the surrounding fluid and generalizes the solution for the hemispheri- 
cal no-slip bump given in an earlier paper by Price (f985). Several expressions including the torque 
and the force acting on the drop will be considered as well as the importance of the deformations 
on the surface for small Capillary numbers. 

PACS numbers: 



INTRODUCTION 



A detailed description of shear flows past viscous droplets resting on a plane wall is necessary for understanding 
and modelling a series of engineering, physiological and natural situations Examples include bubble growth and 
detachment from heated walls @, @f, drop removal and dislocation in detergency problems Q and also membrane 
emulsification, where a liquid is pumped into another, shearing droplets and consequently forming emulsions. Pertinent 
physiological applications can also be included such as the growth of thrombosis and also a tentative description of 
the deformation and dislodging of white and red blood cells adhering to the endothelium. 

Shear flows over protuberances have already been considered in previous works and studies: O'Neill [B| derived an 
infinite series solution for the flow over a full sphere in contact with a wall while Hyman Q and Price Q derived 
a solution for the case of hemispherical no-slip bumps, a problem considered as a limiting case for highly viscous 
droplets. Pozrikidis 0] provided an accurate description of shear flows over protuberances projecting from a plane 
surface usin g th e boundary integral method. A considerable number of papers have also described drop displacement 
problems [9LTT0L [III [I2I [l3| and their variations: Yon, Li and Pozrikidis 14, lj| studied shear induced deformations 
of droplets with fixed contact lines and the work of Schleizer (l6| faced the problem in both pressure driven and shear 
flows. More recently also inertial effects past pinned droplets have been considered 17 1. 

In this paper we derive an exact solution for the case of an incompressible shear flow overpassing a hemispherical 
droplet of arbitrary viscosity. This solution will be computed as a function of the viscosity ratio between the droplet 
and the fluid and is the natural generalization of the one proposed by Price [7j, here recovered in the limit of 
highly viscous droplet. Several analytical expressions including the force and torque acting on the droplet will be 
given. Moreover, in the limit of small Capillary numbers, the role of surface deformations will be investigated and 
comparison with available numerical data will be presented. 

FORMULATION OF THE PROBLEM 

We consider the motion of a viscous incompressible fluid with velocity v, of constant density p and dynamical 
viscosity 77, close to a solid plane boundary. Under the assumption of small Reynolds numbers, this motion is 
reasonably approximated by a uniform shear flow whose magnitude will be denoted with S. We will consider how 
this flow is disturbed by a hemispherical drop of radius R and density p on the plane boundary. We will use polar 
coordinates (r, 9, <p) with the origin in the center of the hemisphere and an x-axis in the direction of the shear flow 
(see figure [T] for details). The flow inside the drop will be identified with a velocity v and viscosity f] in a prescribed 
ratio with respect to the outer viscosity. This ratio will be denoted by A = fj/rj. Indicating with p and p the outer 
and inner pressure respectively, the continuity and Stokes equations read as follows: 

V-v = V-v = 0, (1) 

V-ct = V-<t = 0, (2) 

where 

a = -pl+ T] (Vv + (Vv) T ) , & = -pI + fj (Vv + (Vv) T ) . (3) 
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FIG. 1: Configuration for the problem 



We will assume no-slip boundary conditions on the wall: 



v = r > R, 9 = -, < <f> < 2tt, 



(4) 



v = r<R,9 



< 2tt, 



(5) 



and, on the interface of the hemispherical bump (r = R,0 <9<^fi<(j)< 2n), we should match the inner and 
outer solutions by imposing (i) the kinetic condition for two phases, (ii) the continuity of the tangential velocity in 
two directions and (iii) the continuity of the shearing surface force per unit area in two directions. These boundary 
conditions are written in a general form: 



e r ■ v =e r ■ v = 0, 
e r x v =e r x v, 
e r x {er • e r } =e r x {er • e r } , 



(6) 



with e r the radial normal unit vector at the bump surface. 



GENERAL SOLUTION 



In the subsequent developments we will assume that all variables have been non-dimcnsionalizcd using S, p, rj in the 
outer region and S, p, fj in the inner region together with a reference length scale given by the radius of the drop R. 
We will also introduce the disturbances of the velocity vectors v and v from the linear shear flow as 



q=e x z-\r, q = e x z-\r. 
These disturbances satisfy the continuity and Stokes equations: 

V q = V •(? = (), 



(7) 



(8) 



V 2 q = -Vp V 2 q = -Vp. 



(9) 
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The boundary conditions for non slip walls are now expressed by 



q = r > 1, 6 = ~, < <j> < 2tt, (10) 

q = r < 1, 61 = ~, < < 2tt, (11) 

while, considering the shear velocity e^z written in spherical coordinates, we can rewrite the boundary conditions (J6j) 
using the disturbance velocities as 

q r — cos 8 sin # cos = 0, q r — cos # sin 8 cos = 0, 
90 = qe, qcf = 

dqe , f dq e ~\ dq,p f dq^ A \ (12) 

r = 1, < 6> < tt/2, < < 2vr. 

To reduce the ^-dependence of the boundary conditions (|12p and obtain an axisymmetrical formulation of the problem, 
we introduce u, v, w, it, v and w given by 

q r qe q<t> 
u= -, v = -, w = —— — , 



cos <p cos (p sin q 

q r ~ qe „ 90 
u = -, u = -, w = — 



cos <p cos sin (f> 

We further introduce [/, V and !<¥ as functions of the velocity components u, v and w: 

U =u sin 6 + v cos 8 + w, V = u sin 8 + v cos 8 — w, 
W =u cos 6> — i> sin 8, 

and, for the inner region, U, V and are similarly defined. Then, the boundary conditions on the droplet surface 
can be rewritten using U, V, W, U, V and W as 

Wcos8+ ( U + V ) sinfl-cosflsinfl = 0, (13) 



Wcos8+ [ - + V ) sin 8- cos 8 sin 6> = 0, (14) 



- Wsinfl + £/cos0 = -W' sm8 + Ucos6, (15) 
-Wsin0 + Vcos0 = -Wsin9 + Vcos8, (16) 

^ - 1 ) (-Wsin6» + C/cos6>) = A ( ^- - 1 ) (~Wsm8 + U cos 6>) , (17) 



dr J \dr 



^- - 1 ) (-Fusing + F cos (9) = A -1 ] ( -Wsm8 + V cos (9) . (18) 



^dr J \dr 
As for the pressure field, we define P and P as 



COS (I) COS ffl 
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Since in the creeping flow approximation the pressure is harmonic, we can give a non divergent solution in a series of 
the associated Legendre functions P™: 

00 A P 1 (,A 00 
^=E^T^' P = 5>^V), (20) 

n— 1 n— 1 

where /i = cos# and A n , A n are unknown constants to be determined with the boundary conditions. Here we 
employ the definition of the associated Legendre functions defined in terms of the polynomials of zeroth order as 
P™{x)= (1 - x 2 )™/ 2 d m P„(a;)/da; m (e.g. Pi {co& 6) = sinO). From the Stokes equations (©, we can determine the 
general solutions in the expansion form 

OO A OO f~l 

n=l v ' 11=2 

n=l v ; n=3 



V 



E 2(2n+l)r" P " +lM + E ~ 

-1 V ' n=l 



x 



* = E ( "^2ntir l p "-^^ + E ^P^Cm), 

n=2 ^ ' n=0 

00 4 n+1 00 

* = - E KTIT P «- l(/i) + E ^ n %iW. 

n=3 ^ ' n=l 

* = E "^| 2 ny i p "- l(/i) + E 

n = l ^ ' »=0 

with C„, E n , G n and (7 n , Pn, G n additional constants. Upon using the no-slip boundary condition on the plane wall 
(Pl-CCE]) we get: 

„ _ 2n + l . _ _ (2n + l)A 2n _! 



2(4n+l) zn ' ^ 4(n-l)(4n-l)' 
(2n-l) 2 A 2 „_! . ni 2 „ 

<J2n-l — rj-j TT , <-/2n 



2(4n-l) ' 4n + l 

A _ (n - l)i 2 „-i a _ 2n 2 i 2 „_ 



(2n+l)(4n-l)' 4n - 1 ' 

and, imposing the continuity equation ([5]) together with (fTU|) - (ITT1) . we find: 

_ (2n-l)(2n+l) A 
G2 - 1 " 4(4n-l)(n-l) A2 - l! 

p 1 f (2n-3)(2n+l) 

^ 2n " 2(n-l)(2n-l) \ 4^TI ^ 2 " + &2n 

~ 2n(n - 1) . 

62 - 1 " (2n + l)(4n-l) A2 - 1 ' 



^2n — TT, , , N , n — TT \ : — : — ^2n + Lr2n 



2(n+l)(2n + 3) [ 4n + 1 
while, for the low-order coefficients, we obtain 

A\ = G\ = Pi = Co = G\ = 



and 



A x / 0, Go ^ 0. 

Using all these recursion relations we can rewrite the general solutions for the outer (v 
velocities (v = (v r ,vg,v<p)) representing the solution of the equation set ([I])-©: 



r cos 8 sin 9 — W cos 9 



U + V 



sin ft > cos 



vg = < r cos 2 9 + W sin 9 - 



U + V 



cos ft > cos 



. U V s 
— [ —rcos9 — — + — | sin < 



v r = < r cos 9 sin 9 — W cos t 



sin ft > cos < 



v B = < r cos ft + W sin ft 



tf + v 



COS ft > COS ( 



V4, — I -r cos ft - — + — I sin < 



where 



W = 



10r 2 28r 3 

(n - l)^2n-2 (2n + 1)A 2 » 
(4n - 3)r 2 "- 2 2(4n + l)r 2 ™ 



n=2 



( {2n-l)A 2n -i , (2n+l)(2n + 3)A 2n+1 , x 



V2(4n- l)r 2 



4n(4n + 3)r 2n+1 



ri=2 



28r 3 



2n-2 



2(4n - 3)r 2 "- 2 



(2n-3)(2n+l)A 2n 



G 2 



p2 / 

2(n-l)(2n-l)(4n+l)r 2 ™ ' 2{n - l)(2n - l)r 2n ' " 1 



/ A 2n _! (2n + 3)A 2 » +1 \ 2 

^2(4n - l)^"- 1 4n(4n + 3)r 2 "+! / 2n W 



r 



9A 3 P 2 (//) 

+ 2_, [{ (4n - 3)r 2 - 2 + T 2 ^ J P2 - l(Ai) 

n(2n - l)A 2n -i _ (2n + l) 2 A 2ra+ i \ . 
' (4 n _i) r 2„-i 2(4n + 3)r 2 "+! J 2 "^ j 
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W 



3A 2 r 3 Pl{u) 
10 

^( (2n+l)A 2n r^ (n - l)A 2n . 2 r^ \ ,„ 



n=2 



2(4ra+l) 



4n - 3 



A / (n + l)i 2n+1 r 2 "+ 2 2n(n - l)i 2 „_ir 2 " , pl . . 
+ ^ i o h TnZ. , TT77I 7T~ I ^2nW> 



n=l 



4?i + 3 



(2n + l)(4n- 1) 



(30) 



i 2 „r 2 " +1 _ 2(n-l)(n + l)i 2n _ 2 r 2 "- 1 
2(4n+l) n(2n + l)(4n-3) 



G 2n - 2 r 2n ~ 
2n(2n + 1 

y ( i 2n+ ir 2 "+ 2 (n - l)i 2 „-ir 2 " 

71=1 V 



^2n-l(M) 



2(4n + 3) (2n + l)(4n- 1) J P * n ^> 



(31) 



s Z^ 



n(2n+ l)i 2 „r 2 " +1 
4n + 1 



+ G 2 „_ 2 r 2 "- 1 P an _ x 0i) 



(n+l)(2n+l)i 2n+1 r 2 "+ 2 2n 2 i 2 „_ 1 r 2 ™ 



4n + 3 



An - I 



(32) 



In the previous expressions the six unknown sets of coefficients A 2n+ \, A 2n , G 2n , A 2n —i, A 2n and G 2n - 2 with n > 1 
have now to be determined by the imposition of the six conditions (|13p - (|18p on the hemispherical bump. 



Determination of the coefficients 

To determine the coefficients A 2n+ \, A 2n , G 2n , A 2n _\, A 2n and G 2rl _ 2 for n > 1, we need recurrence formulae that 
can be deduced from the boundary conditions at r = 1 given by (fT3|) - (fT8|) . These recurrence formulae are given in the 
appendix and, although they are written as infinite series, we numerically solve the problem by truncating the system 
of equations determining the coefficients. For each of the six coefficients we consider modes up to n = 100 and to 
have a good indication of the accuracy of this procedure, the numerically estimated coefficients are substituted into 
(Ti"3")) - (|18p until the boundary conditions are satisfied to at least five decimal places. The coefficients A 2n+ i, A 2n , G 2n , 
A 2n -i, A 2n and G 2n - 2 up to n — 15 are listed in tables HWl for A = 0, 0.1, 1, 10 and 10 100 and to obtain the full 
solutions, they must be used into equations ([2"T j) - (f32"|) [2^ ]. 

In figure [2] we show the kinematical structure of the flow and the development of stagnation points and flow reversal. 
Let us notice that for high viscosity ratio the streamlines plot is consistent with what obtained by Pozrikidis 
It is now useful to check the validity for two extreme limits A — > and A — > oo. For the free-slip bump, A = 0, as also 
confirmed by our numerical solution shown in Table HI the solution is analytical and can be expressed as 

2 2 

A n = — S n2 , G n — —S n2 . 

3 5 

Using this solution, we determine U, V, W and P 



The velocity field around the hemisphere is then given by 



v =e, r I r I cos 8 sin 9 cos (f> + egr cos 2 9 cos (j> — e^r cos 9 sin < 



-2-1012-2-1012 




FIG. 2: Streamlines patterns in the in xz mid-plane (for y = 0). 



that can also be rewritten together with the pressure in Cartesian coordinates 

(33) 



xyz 

v =e x I z — I - e y — = e z — =- 



2xz 

It should be noted that the present solution Q33p is same as the flow around a free-slip sphere in an infinite fluid. It 
is because the presence of the free-slip sphere does not alter the 9 and (f> components of the velocity vector from the 
uniform shear flow and then the modulated velocity becomes zero at z = 0. 

In the case of the largest viscosity ratio A = 10 100 , the coefficients A2 n +i, and Gn n listed in Table fVl show 
quantitative agreement with those for the no-slip bump shown in the paper by Price [7|. The inner solution is 
analytical and given by A n = and G n = 25o n that imply a zero velocity inside the bump. 



Torque and Force 



We now evaluate the force and the torque acting on the hemisphere for various viscosity ratios. The general strategy 
will be to write them as a functions of the constants A2 n , S^n+i, Gin and then use the coefficients obtained from the 
numerical solution in the previous section in order to evaluate them. To do this, we write the force F and the torque 
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T in the most general form 



F = I d<j> d9 sm9 (<x ■ e r ) r=1 , 
Jo Jo 

T = I dcf) I dO siri9 [re r x (cr • e r )] r= i. 
Jo Jo 



Let us notice that the symmetry of the system with respect to 4> is such that only the x component of the planar force 
vector (F x ,F y ) and y component of the torque vector T y are non-zero. To estimate these components, considering 
the kinetic condition in (|12[) and using U, V and W given by the general solution we can exactly express the torque 
and force as a function of A2 n ,Am+\ and G2n- 



Fr 



E 



(-l)"(2n + l)!! 



5" u ^ (2n-3)(2n-l)(2n-2)!! 
(4n 3 - 9n - 2)A 2n (2?i 2 - 3n - l)G 2n 



2(n + l)(4n + l) 



n(2n+ 1) 



(34) 



with 



Ty =1T 



3 A 2 + G 2 + ± ^^ A 2 ^ 



^ (2n-l)(2n)! 



(35) 



(2n + 1)!! = (2n + 1) x (2n - 1) x ... x 1, 



(2n)!! = (2n) x (2n - 2) x ... x 2, 0!! = 1. 

In the limit of highly viscous droplet (A — > oo), from the numerical evaluation of the coefficients for A = 10 100 (see 
Table fV]) we determine the force and torque as 

F x = 4.30322tt, T y = 2.44132tt, 

and both results are consistent with the evaluations F x = 4.30^,^ = 2.44-7T for the rigid no-slip bump given in the 
paper by Price Q- 

The profiles of F x and T y as a function of A are respectively shown in figure [3] and figure [U As is clearly shown both 
the torque and the force approach a finite value for very small (A — > 0) and large (A — * oo) viscosity ratio. Considering 
these asymptotic behaviors, we can estimate 



F x _2 + 4.51003A 
~ ~1 + 1.04806A' 



T 



2.18808A 



7T 1 + 0.896271A' 

in which the numerical coefficients have been determined from the method of least squares. 



(36) 
(37) 



Effects of Capillary number on deformation 

In this section we study the local difference of the normal stress between the inner and outer regions of the bump. 
If we assume a finite surface tension 7 at the droplet surface, the difference between the the two normal stresses 
must equal the surface tension times the local curvature of the hemispherical bump. In other words, for fixed surface 
tension, this gives the possibility to estimate the deformations on the droplet surface consistently with the Laplace 
law [II Ql]. Resting on the assumption that these deformations are small, we can introduce a Capillary number as 



Ca 



r/SRo 



9 



4.5 
4.0 
3.5 
^ 3.0 
2.5 



2.0 o ooooooooo -e 




e oooooooooo 



10" 4 10" 3 10" 2 10" 1 10° 10 1 10 2 10 3 10 4 

A 

FIG. 3: The force acting on the hemisphere versus the viscosity ratio. The symbol o corresponds to the evaluation based on 
(|34|) . The solid line corresponds to the approximation (|36|l . 



where Rq, S and 7 are dimensional length, shear and surface tension scales. Next we write a dimensional stress tensor 
as 



(dim) 



crS-q - P s , 



(dim) 



&Si]X - P s 



where the contributions er and <r are purely dimensionless, i.e. they correspond to a unitary shear S and unitary 
viscosity r\ described in the previous sections. The terms P s and P s represent constant dimensional background 
pressures, i.e. the pressure inside/outside the bump in absence of flow. Assuming Ca <C 1, if we indicate with r = R 
the position of the interface, we should write: 



R{0, <j>) = Rq(1 + CaR^{9, <£) + 0(Ca 2 )), 



(38) 



where R^ accounts for the displacement from the pure hemisphere. The Laplace law on the interface is written in a 
general form 

K (dim)^ _ a (dim) _ fi(dim) ^ 



where K^ dlm ) = Ji- denotes the dimensional curvature, with k its dimensionless counterpart. If Ca <C 1, the displace- 
ment is small enough and the curvature on the bump interface is given by 



k =R < V 



V(r - R) 
Mr-R)\ 



R {V 2 (r-R)} 



r=R 



l + CoR« - 
-2Ca° ~ Ca 1 



1 d ( . 8 
smfci— 



sin 9 86 



r=R 

1 d 2 



1 d ( . d_ 

sin 86 [ 00 



sm c 

1 d 2 



CaR {1 



sm 



+ 0(Ca 2 ), 
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10" 4 10" 3 10" 2 10" 1 10° 10 1 10 2 10 3 10 4 



FIG. 4: The torque acting on the hemisphere versus the viscosity ratio. The symbol o corresponds to the evaluation based on 
(|35|) . The solid line corresponds to the approximation (I37|l . 



hence ^ 



2 

Chi 



(P s - P.)/(Srj) and 
1 d 



sin 6 80 V 30 



1 d 2 



sin 2 6 < 



(39) 



: (a rr - Aay r ) r=1 + (P s - P s )/(Sr)). 



In the absence of the shear flow, the pressure is constant inside and outside the bump and equal to P s and P s 
respectively. The gap between the two pressures is determined by (P s — P s )/(Srj) = 2Ca _1 from the C(Ca _1 ) of 
(|3"9")) . This means that in presence of a non zero velocity, in order to find consistency with the Laplace law, we should 
consider a finite displacement given by a non zero R^ in (|38[) . To determine this displacement, equation (|3"9")l is 
solved at O(Ca ) : 



with the boundary condition 



1 d ( . 8 



= at 



1 



if 



sin 2 9 > 



R m = (a r 



, < 6 < 2?r. 



(40) 



(41) 



This condition corresponds to fix the contact line while accounting for the shape deformation of the droplet. More 
generally, the contact line will be deformed and in principle one can solve for the contact line shape as part of the 
problem [12|)[l3j]- However, the requirement ((4T|) is still useful, given the comparison it allows with other works [HI, EH]. 
If we write the normal stress jump in the expansion form with respect to the associated Legendre polynomials 



A<7 rr .) r=1 = ^ Q n Pn(n) COS( 
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the expansion coefficients Q n can be determined as: 

Q 2n =^(l-A) + 4n2 + 6n - 1 A 2 „ 
V2n 3 v ' 4n-l 



2(n + l) 2 (2n + 3)(8n 2 + 6n - 3) 2(n + 1) 

+ n(2n+l)(4ra + 3)(4n + 5) 2 " +2 + n(2n + 1) 2 " +2 

/ 2(n - l)(2n - l) 2 (4n 2 + n - 2) > (42) 
+ ^ n(2n+l)(4n-3)(4n- 1) 2 "~ 2 
4n 2 - 2n - 3 - 2n - 1 - \ 

4n + 3 n(2n + l) / 



4n 2 + 2n-3 (n + l)(2n+l) 2 

An - 3 n(4n + 1) 

^ 4(n-l) 2 (2n-3) > 4n 2 - 6n-l , ' 

+ A ( v (2n-l)(4n-3) Mn ~ 3 + 4n + l 



We now solve (|40l) and determine the displacement 



n=2 

where 



flW = p + £ fl n )p n(M) ] COS0, (44) 



ff (D y (-l)»(2fc-l)ll 

^ -Z. 2 (fc_l)(2A ; + l)(2fc-2)!! g2fc - 1 ' (45) 



k=2 



^- ( »-ix» + » for " a2 - (46) 

Concerning the dependence of the normal stress jump, this is simply proportional to cos0. In figure [5] we show the 
deformation of the bump on the plane of y — with a Capillary number Ca = 0.05. These deformations have been 
estimated from fH]) . (|45[) an d PS)) with the numerical coefficients in (|12"1) and P3"| determined with the procedure of 
the previous section. Moreover, to better quantify the role of deformations with respect to the hemispherical bump, 
advancing and receding contact angles should be respectively defined by w/2 + A9 and ir/2 — A6, with Ad proportional 
to Ca if the deformation is small enough (see figure [H]). 

For various viscosity ratios A, A9/Ca is then plotted in figure [7l Considering the asymptotic behaviors, this 
deformation can be approximated by 

AO « (2 + 7.07185A) 

(1 + 0.908826A) ' v ' 

in which the numerical coefficients have been determined from the method of least squares. 

For the case of hemispherical droplet whose viscosity is equal to the surrounding fluid (A = 1), it is instructive to 
compare our results with the ones of Li and Pozrikidis [3] ■ 

In figure [8] we plot the deformation parameter D = ^+3 as a function of the Capillary number. A and B are, 
respectively, the maximum and minimum radial distances of the interface from the ori gin . Note that the corresponding 
slope for a spherical drop suspended in an infinite simple shear flow is equal to 35/32 [2fJ and suggests that the wall is 
promoting the deformation of the drop as it is also confirmed by the computation of Li and Pozrikidis [lij ] (see their 
figure 6 (e)) with the boundary integral method. To further compare our results with the ones of Li and Pozrikidis, 
in figure [9] we plot the advancing and receding contact angles as a function of the Capillary number and a good 
agreement is found. 

The numerical investigations of Li and Pozrikidis (lij were also extended in another paper by Yon and Pozrikidis [IH 
by examining the role played by the viscosity ratio. In particular in figureHUlwe compare our result for the advancing 
contact angle as a function of the Capillary number with the numerical results shown in figure 3 (g) of the paper by 
Yon and Pozrikidis and quantitative agreement is found. 
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FIG. 5: Shape of the bump on the plane of y = for various viscosity ratios A. The dashed line shows the hemisphere with 
no deformation. The solid line shows the deformed interface for a Capillary number Ca — 0.05. 




FIG. 6: Advancing and receding contact angles due to deformations are introduced in the problem in the small Capillary 
number regime. 



CONCLUDING REMARKS 

An exact solution for the problem of a linear shear flow overpassing a hemispherical droplet resting on a plane wall 
has been explicitly computed. The results here presented are the natural generalization of the ones proposed in an 
earlier paper by Price [7[ who investigated the case of a hemispherical no-slip bump. Several expressions including 
the torque and force acting on the hemisphere are exactly derived and computed as a function of the viscosity ratio. 
The small deformations on the hemispherical surface have also been quantified in the limit of small Capillary numbers 
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10" 4 10" 3 10" 2 10" 1 10° 10 1 10 2 10 3 10 4 



FIG. 7: Contact angle deformations as a function of the viscosity ratio in the limit of small Capillary numbers. The symbol o 
corresponds to the solution based on (|40|l . The solid line corresponds to the approximation (|47l) . 




FIG. 8: Deformation parameter D as a function of the Capillary number for a viscosity ratio A = 1. Our solution (straight line) 
is compared with the results obtained with the boundary integral method (o) by Li and Pozrikidis (1996). The deformation 
parameter for the case of a droplet suspended in an infinite shear flow is also reported (dashed line). 
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FIG. 10: Advancing contact angles as a function of the Capillary number for viscosity ratios A = 0, 1, 10. Our results are 
compared with the numerical results obtained with the boundary integral method (o) by Yon and Pozrikidis (1999). 
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and a comparison with available results has been presented. Extensions of the present study can be made into several 
directions. One extension concerns the generalization of the present results to non hemispherical droplets [§], Il5j |. 
In this case, one would apply the same procedure discussed in the paper but the equation set emerging from the 
imposition of the no-slip boundary on the plane wall needs to be solved numerically. Another interesting extension 
would be the possibility to study the transient deformation and asymptotic shapes of drop induced by insoluble 
surfactant and also the effect of gravity. This is another case where it would be important to complement other 
numerical works present in the literature 13, 13, 14, T^. Finally, all results presented in this paper are also useful in 



the context of ensemble averaging techniques to formulate effective boundary condition problems 21, 23, 23. 24. 2 
an issue which has been recently renewed due to its relevance in small scale hydrodynamics [2~ 



APPENDIX 



In this appendix we report the recurrence formulae to determine the coefficients A2 n +i, Am-, Gim Azn—i, A^n and 
Gi n -i for n > 1. The relevant boundary conditions at r = 1 are given by l]13 p -(|18 |) . Using the general representations 
(|27[) - (|32[) . we can rewrite the boundary conditions and apply definite integrals for Legendre polynomials 



1 

dfl P2n(p)P2k-l{^) 



Jo 



1 

2 



6 nk (2fc-l + m)! 
(4fc-l) (2fc-l-m)!' 

S nk (2k + m)\ 
(4k + 1) (2fc-m)!' 

(_l)"+fc(2n-l)!!(2fc-l)!! 
~ 2(2n - 2k + l)(n + fc)(2ra)!!(2fc - 2)!! ' 

(_l)"+fc(2n + l)!!(2fc- 1)!! 
2(2n -2k + l)(n + k)(2n - 2)!!(2fc - 2)!! ' 

(_l)"+fc(2n + l)!!(2fc + l)!! 



/„ ' 2(2n-2fc+l)(n + fc)(27i-2)!!(2fc-4)!!' 

In this way, we can eliminate the summation for the even-order mode and we finally obtain the recurrence formulae: 

hi _ (2fc + 1) (k + l)(2fc + 3)(8fc 2 + 6k - 3) 

T ~2(4fc- 1) 2k + 2k(2k + l)(4fc + 3)(4A: + 5) 2k+2 

+ 2k(2k + l) G2k+2 ~ Afc ^^ Tl + l )^*Mn+\ 

^ ' n— 1 

S kl (fc - l)(2fc - l)(4fc 2 + k-2) ■ k - 

3 ~ fc(2fc + l)(4fc-3)(4fc-l) 2k ~ 2 + 4k + 3 2k 

1 X 

+ 2fc(2fc + l) ^ fc - 2 - Ak 5>» - ^M^, 

x ' n— 1 
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6fc 2 -7k -2 (2fc + 3)(6fc 2 + 5fc-3) 

-A2fe + o; , o; — ; — ; ,,,,, — — rr^2fe+2 



(2fc-l)(4fc-l)(4fc + l) 2fc(2fc + l)(4fc + 3)(4fc + 5) " 

I 1 r | (2fc + 3) 

+ (2fc-l)(4fc-l) 2fe 2fc(2fc + l)(4fc + 3) 2k+2 

A °° 

$^(n - l)(2n + 3)A„, fcJ 4 2n+1 



(2fc-l)(2fc + 2) 



(fc - l)(6fc 2 + /c - 4) j 6fc 2 + 13fc + 3 

-^2fe-2 — 7777 — , 77777 — , 77771 — 7~r77^2/c 



fc(2fc + l)(4fc-3)(4fc- 1) 2(fc + l)(4fc + l)(4fc + 3)' 

Glk-2 + 7771 — , 7T771 — — 777^*2*; 



fc(2fc + l)(4fc- 1) 2(fc + l)(4fc + 3) 

A 00 



2fc 2 (2fc + l) (fc + l)(2fc + l)(2fc + 3) 

r^2fe ... ; w,, — — TT ^2fe+2 



(4fc- l)(4Jfe + l) (4fc + 3)(4fc + 5) 

2fc 2k + 1 

+ 4fc - 1 G 2 fc + 4fc - 3 G 2 fc +2 + A fe 22 n ( 2n + l)K t ,kA 2n +i 

2fc(fc- l)(2fc- 1) - fc(2fc+l) 2 ~ 

1 (4fc - 3)(4fc - 1) 2fc ~ 2 + (4fc+l)(4fc + 3) 2k 
2k 2k + 1 °° 

n=l 



(2fc + l)(6fc 2 -7fc-2) (2fc + 3) 2 (6fc 2 + 5fc-3) , 

7^2fc - n , , n , ; : w,, ; w,, , r x ^2fc+2 



(2k- l)(4fc- l)(4fc + 1) 2fc(2fc + l)(4fc + 3)(4fc + 5) 

(2fc + l) (2fc + 3) 2 

-<-f2fe — 7777771 — 77771 — r _ r77 ( - T 2fc+2 



(2fc-l)(4fc-l) 2fc(2fc + l)(4fc + 3) 

A 00 

+ (fc + i)^-!) 5> - + l ^ 2n + 3 ) A »,^+i 

f 2(fc- l) 2 (6fc 2 + fc-4) ~ fc(6fc 2 + 13fc + 3) , 

| yfc(2fc + 1)(4A; - 3)(4fc - 1) 2fe ~ 2 ~ (k + l)(4fe + l)(4fc + 3) 2fc 

2(fc-l) 2 > fc > 1 

fc(2fc + l)(4fc - 1) 2fc_2 (fc + l)(4fc + 3) 2k [ 



+ (fc + ^ _ 1} X> - X )( 2n - !)( 2n + 3 ) A n,fc^ n . i ■ 
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where we have used 



and 



2fc 2 (2/c+l) 2 (fc+l)(2/c + l)(2fc + 3) 2 

-^2fc H ,., ; rrr— — — rr ^2fe+2 



(4fc- l)(4fc + l) (4/c + 3)(4fc + 5) 

2fc(2fc + l) (2fc + l)(2fc + 3)^ 

" (Ak-l) ° 2k (4fcT3) &2k+2 

OO 

- 2A k n(n + l)(2n + l)A„^A 2 „+i 

n=l 

J 4fc(fc — l) 2 (2fc — 1) j 2fc 2 (2fc+l) 2 ~ 

;A j (4fc - 3)(4A - 1) 2fe " 2 + (4TTT)(4fcT3) 2fe 

4fc(fc-l) ^ 2fc(2fc+l) ^ 1 

+ 4fc-l G2fe - 2 + 4fc + 3 G2fe j 

OG 

- 2AA fc ^ n(2n - l)(2n + l)A nife i 2 „_i, 

n=l 



(-l)"(2n + l)!! 



A 



2(2n - 2fc - l)(2n - 2fc + l)(n + k)(n + k + l)(2n)!! ' 

(-l)"(2n-l)!! 

~4(2n -2k- l)(2n - 2fc + l)(n + fe)(n + fc + l)(2n - 2)!! : 

_ (-l) fc (4fc + l)(2fc- 1)!! 
k ~ (2kJ\ ' 
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TABLE I: Coefficients for A = 0. 



n 




^2n + l 






Mn 








G2n 








A-2„ 




1 


-4. 


.2435 


X 


10" 


12 






0.66667 






0.40000 


— 11.47649 


2.53882 


8.03099 


2 


-2 


.3846 


X 


10" 


■12 


-9 


.0339 


X 


io- 


■14 


-2 


.5928 


X 


io- 


13 


-2.13519 


0.17201 


5.42824 


3 


1 


.1579 


X 


10" 


■12 


-1 


.8727 


X 


io- 


■14 


2 


.8944 


X 


10" 


13 


0.24983 


-0.071796 


0.12475 


4 


— 1 


.8929 


X 


io- 


■12 


-7 


.2444 


X 


io- 


■14 


— 1 


.8123 


X 


io- 


■13 


—0.062825 


0.036175 


-0.082789 


5 


—2, 


.1743 


X 


10" 


■13 


1 


.6780 


X 


io- 


■14 


2 


.5391 


X 


io- 


■13 


0.019726 


-0.020377 


0.056309 


6 


-2 


.7782 


X 


10" 


■13 


— 1 


.5622 


X 


io- 


■14 


-2 


.1449 


X 


io- 


■13 


-0.0063530 


0.012379 


-0.039535 


7 


3 


.3070 


X 


io- 


■15 


1. 


.3030 


x 


10" 


■14 


1 


.6910 


x 


10" 


■13 


0.0015681 


—0.0079442 


0.028544 


8 


—8. 


.3398 


x 


10" 


■14 


4 


3793 


X 


io- 


■15 


—2 


.0627 


X 


io- 


■13 


0.00024267 


0.0053167 


-0.021091 


Q 

y 


7 


.9357 


X 


io- 


■15 


o 
o 


.6435 


X 


io- 


■14 


y 


.0663 


X 


io- 


■14 


— u.uuuyuoyy 


— U.UUOO ( 


U.UlOOOi 


10 


-1. 


.8344 


X 


io- 


■13 


-1 


.1164 


X 


io- 


■14 


-1 


.4936 


X 


io- 


■13 


0.0011024 


0.0026142 


-0.012147 


11 


1 


.7744 


X 


10" 


■14 


1 


.6463 


X 


10" 


■14 


1 


.3542 


X 


10" 


■13 


-0.0011030 


-0.0018995 


0.0094079 


12 


-3. 


.0430 


X 


10" 


■13 


6 


.7278 


X 


10" 


■15 


-2 


.5442 


X 


10" 


■13 


0.0010260 


0.0014061 


-0.0073615 


13 


-2. 


.2858 


X 


10" 


■13 


1 


.5070 


X 


10" 


■14 


7 


.3790 


X 


10" 


■14 


-0.00092290 


-0.0010573 


0.0058075 


14 


-2. 


.6355 


X 


io- 


■13 


-1 


.8224 


X 


io- 


■14 


-5 


.9443 


X 


io- 


■14 


0.00081663 


0.00080565 


-0.0046107 


15 


2 


.3231 


X 


io- 


■13 


-1 


.0819 


X 


10- 


■14 


2 


.0463 


X 


io- 


■13 


-0.00071683 


-0.00062094 


0.0036777 



TABLE II: Coefficients for A = 0.1. 



n 


A2n + 1 


A 2n 


G2n 


A.2n-1 


A 2n 


G2n-2 


1 


-0.17436 


0.84972 


0.50983 


-10.46143 


2.31835 


7.49147 


2 


-0.077633 


0.15800 


-0.077149 


-1.94082 


0.15400 


4.92646 


3 


0.012187 


0.0094929 


-0.0011555 


0.22749 


-0.064461 


0.11746 


4 


-0.0035702 


-0.0046015 


0.0023703 


-0.057379 


0.032545 


-0.077592 


5 


0.0012341 


0.0025121 


-0.0022646 


0.018100 


-0.018362 


0.052686 


6 


-0.00042666 


-0.0014877 


0.0019020 


-0.0058827 


0.011167 


-0.036976 


7 


0.00011347 


0.00093504 


-0.0015382 


0.0014964 


-0.0071731 


0.026705 


8 


0.000013522 


-0.00061510 


0.0012302 


0.00017241 


0.0048045 


-0.019748 


9 


-0.000063680 


0.00041948 


-0.00098250 


-0.00079079 


-0.0033262 


0.014888 


10 


0.000080443 


-0.00029453 


0.00078655 


0.00097802 


0.0023654 


-0.011407 


11 


-0.000082443 


0.00021183 


-0.00063199 


-0.00098504 


-0.0017197 


0.0088482 


12 


0.000078085 


-0.00015543 


0.00050977 


0.00091948 


0.0012736 


-0.0069382 


13 


-0.000071276 


0.00011601 


-0.00041261 


-0.00082899 


-0.00095812 


0.0054870 


14 


0.000063845 


-0.000087837 


0.00033493 


0.00073472 


0.00073043 


-0.0043687 


15 


-0.000056630 


0.000067335 


-0.00027242 


-0.00064571 


-0.00056322 


0.0034961 
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TABLE III: Coefficients for A = 1. 



n 




A2n 


Cr2n 


A^n 1 




C*2n 2 


1 


—0.97542 


1.68344 


1.01007 


-5.85249 


1.31294 


5.05033 


2 


-0.42739 


0.87935 


-0.40999 


-1.06849 


0.077130 


2.67984 


3 


0.067950 


0.048954 


-0.014276 


0.12684 


-0.033073 


0.077672 


4 


-0.020266 


-0.024073 


0.017863 


-0.032570 


0.016955 


-0.049837 


5 


0.0071884 


0.013276 


-0.015674 


0.010543 


—0.0096669 


0.033457 


6 


—0.0026054 


—0.0079212 


0.012751 


-0.0035922 


0.0059251 


-0.023387 


7 


0.00079580 


0.0050080 


-0.010168 


0.0010494 


-0.0038295 


0.016889 


8 


-0.000042750 


—0.0033104 


0.0080860 


-0.000054498 


0.0025779 


—0.012520 


Q 
y 


— U.UUUZDojZ 


u.uuzzooy 




n nnn^^^99 

— U.UUUOOOZZ 


— u.uui ( yzo 




10 


0.00038388 


-0.0015974 


0.0051775 


0.00046672 


0.0012795 


-0.0073022 


11 


-0.00041078 


0.0011527 


-0.0041798 


-0.00049081 


-0.00093345 


0.0057059 


12 


0.00039788 


-0.00084839 


0.0033941 


0.00046853 


0.00069351 


-0.0045116 


13 


-0.00036830 


0.00063500 


-0.0027709 


-0.00042837 


-0.00052328 


0.0036025 


14 


0.00033312 


-0.00048214 


0.0022727 


0.00038336 


0.00040007 


-0.0029001 


15 


-0.00029762 


0.00037060 


-0.0018715 


-0.00033936 


-0.00030935 


0.0023503 



TABLE IV: Coefficients for A = 10. 



n 


A-2n+l 


A 2n 


G2n 


A2n-1 


A 2 n 


G2n-2 


1 


-1.83516 


2.55162 


1.53097 


-1.10110 


0.25418 


2.56548 


2 


-0.78449 


1.64184 


-0.68804 


-0.19612 


0.011796 


0.47557 


3 


0.12819 


0.079556 


-0.050045 


0.023929 


-0.0054338 


0.018141 


4 


-0.039327 


-0.040886 


0.045731 


-0.0063205 


0.0028862 


-0.011069 


5 


0.014454 


0.023098 


-0.036923 


0.0021199 


-0.0016810 


0.0072831 


6 


-0.0055521 


-0.013997 


0.029024 


-0.00076551 


0.0010455 


-0.0050478 


7 


0.0019496 


0.0089486 


-0.022775 


0.00025708 


-0.00068302 


0.0036351 


8 


-0.00039965 


-0.0059665 


0.017979 


-0.000050954 


0.00046375 


-0.0026961 


9 


-0.00027530 


0.0041147 


-0.014313 


-0.000034189 


-0.00032475 


0.0020468 


10 


0.00055414 


-0.0029172 


0.011495 


0.000067373 


0.00023324 


-0.0015835 


11 


-0.00064763 


0.0021163 


-0.0093099 


-0.000077380 


-0.00017109 


0.0012443 


12 


0.00065388 


-0.0015652 


0.0075980 


0.000076998 


0.00012775 


-0.00099055 


13 


-0.00062032 


0.0011769 


-0.0062438 


-0.000072149 


-0.000096843 


0.00079718 


14 


0.00057041 


-0.00089748 


0.0051623 


0.000065643 


0.000074373 


-0.00064750 


15 


-0.00051578 


0.00069277 


-0.0042910 


-0.000058812 


-0.000057759 


0.00053005 



TABLE V: Coefficients for A = 10 



n 


A2n + 1 


A 2n 


G 2 „ 


A 2 


n- 


-1 






A 2n 






G 2 


n - 


-2 




1 


-2.04530 


2.75543 


1.65326 


-1.2272 > 


( 10 


-99 


2 


.8662 


X 


10" 


100 






2.00000 


2 


-0.86783 


1.82575 


-0.73104 


-2.1696 


X 


io- 


100 


1 


.2213 


X 


10" 


■101 


5. 


1794 


X 


10" 


■100 


3 


0.14321 


0.084338 


-0.063687 


2.6733 


X 


io- 


101 


-5 


.8109 


X 


10" 


-102 


2. 


1417 


X 


10" 


-101 


4 


-0.044300 


-0.044198 


0.054743 


-7.1196 


X 


io- 


102 


3 


.1280 


X 


10" 


-102 


-1. 


.2864 


X 


10" 


-101 


5 


0.016436 


0.025196 


-0.043364 


2.4106 


X 


io- 


102 


-1 


.8352 


X 


10" 


-102 


8. 


.4124 


X 


10" 


-102 


6 


-0.0064047 


-0.015350 


0.033811 


-8.8306 


X 


io- 


103 


1. 


.1467 


X 


10" 


-102 


-5. 


.8150 


X 


10" 


-102 


7 


0.0023185 


0.0098491 


-0.026427 


3.0574 


X 


io- 


103 


-7 


.5171 


X 


10" 


-103 


4. 


1832 


X 


10" 


-102 


8 


-0.00054568 


-0.0065846 


0.020824 


-6.9572 


X 


io- 


104 


5 


.1171 


X 


10" 


-103 


-3 


.1021 


X 


10" 


-102 


9 


-0.00023620 


0.0045508 


-0.016566 


-2.9333 


X 


io- 


104 


-3 


.5909 


X 


10" 


-103 


2. 


.3560 


X 


10" 


-102 


10 


0.00056676 


-0.0032322 


0.013305 


6.8908 


X 


io- 


104 


2 


.5837 


X 


10" 


-103 


-1. 


.8242 


X 


10" 


-102 


11 


-0.00068450 


0.0023485 


-0.010782 


-8.1786 


X 


io- 


104 


-1 


.8982 


X 


10" 


-103 


1. 


.4350 


X 


10" 


-102 


12 


0.00070095 


-0.0017395 


0.0088074 


8.2541 


X 


io- 


104 


1 


.4194 


X 


10" 


-103 


-1. 


.1439 


X 


10" 


-102 


13 


-0.00067033 


0.0013096 


-0.0072468 


-7.7965 


X 


io- 


104 


-1 


.0775 


X 


10" 


-103 


9 


.2200 


X 


10" 


-103 


14 


0.00061962 


-0.00099993 


0.0060008 


7.1306 


X 


io- 


104 


8 


.2854 


X 


10" 


-104 


-7. 


5021 


X 


10" 


-103 


15 


-0.00056237 


0.00077279 


-0.0049970 


-6.4124 


X 


io- 


104 


-6 


.4425 


X 


10" 


-104 


6. 


.1534 


X 


10" 


-103 
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